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^ \ Abstract 

■ Subtle problems in the theory of neutrino oscillations in vacuum are discussed Q. 
It is shown that Lorentz invariance implies that in general flavor neutrinos in oscillation 

\ experiments are superpositions of massive neutrinos with different energies and different 

^ ' momenta. It is argued that a wave packet description of massive neutrinos is necessary 

■ in order to understand the physics of neutrino oscillations. 
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0\ , 

^ 1 Introduction 

m ; 

Q \ The standard theory of neutrino oscillations in vacuum p[ has been developed in the 70's 
^ I on the basis of four main assumptions: 

^ , (Al) Neutrinos are extremely relativistic particles. 

^ , (A2) Neutrinos produced in charged-current weak interaction processes together with 

(D ' charged leptons a+ are described by the flavor state 

^ ■ V- 

>: W) = 2_^Kt,\vk) (a = e,/i,r), (1) 

• 1— I , I. 
>< ■ 

^ • where U is the mixing matrix of the neutrino fields, v^l = ^k^ak^kL- The massive 

■ ■ neutrino states \h'k) are quanta of the fields z/^ of neutrinos with mass m^. Their energy 

Ek and momentum pk are connected by the relativistic dispersion relation 

El=pl + ml. (2) 

Since in oscillation experiments neutrinos propagate along a macroscopic distance be- 
tween production and detection, we consider only one spatial direction along the neu- 
trino path. 

(A3) The massive neutrino states have the same momentum p ("equal momentum as- 
sumption", pk =p), but different energies. 



2 2 

■mt „ mi 



Ek = ^P^ + ml^p+^ = E + ^. (3) 

Here E = p the energy of a massless neutrino and the approximation is valid for 
extremely relativistic neutrinos. 
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(A4) The time T of propagation of neutrinos from source to detector is approximately equal 
to the source-detector distance L. 

The massive neutrino states |z/fc) evolve in time according to the Schrodinger equation, whose 
solution is 



Hit)) 



-iEi,t 



(4) 



Using the assumptions |(A1)| - [[A4)1 it is straightforward to derive the oscillation probability 



Pua^U/s (T) 

as a function of the time T and the oscillation probability 



k 



(5) 



P..^.„iL) = \U^k\'\U^k\' + 2Re KkUpkU^jU;^ exp f -.^'"'^•'■^ 



k 



k>j 



2E 



(6) 



as a function of the distance L, with Am|^. = m\-~rn?^. The probability (H) is the one observed 
in real experiments where the propagation time T in not measured but the source-detector 
distance L is known. 

Let us examine critically the assumptions |(A1)|-[(A¥). 



The assumption KA1)| is easily understood to be correct by noticing that, although not 
all existing neutrinos are extremely relativistic, in neutrino oscillation experiments one is 
interested in detectable neutrinos, which have energy larger than some fraction of MeV. 
Indeed, neutrinos are detected in: 

1. Charged current weak processes which have an energy threshold larger than some 
fraction of MeV. For example^: 

— Eth = 0.233 MeV for + '^Ga ^^Ge + e' in the GALLEX [|, SAGE [| and 
GNO P solar neutrino experiments. 

— i?th = 0.81 MeV for v^, + ^''Cl — >■ ^''Ar -|- e~ in the Homestake solar neutrino 
experiment. 

— i?th = 1-8 MeV for z/g + p — n + e"*" in reactor neutrino experiments (for example 
Bugey and CHOOZ [|). 

— Eth = 110 MeV ioi p^ + n ^ p + fi'. 

— Eth — m^/2me = 10.9 GeV for z/^ + e~ ^ Ue + fi~ . 



2. The elastic scattering process u + e" — > u + , whose cross section is proportional 
to the neutrino energy {(t{E) ~ aQE/rrif., with ctq ~ lO'^^^cm^). Therefore, an energy 
threshold of some MeV's is needed in order to have a signal above the background. 
For example, Eth — 5 MeV in the Super-Kamiokande solar neutrino experiment. 



^ In a scattering process v + A ^ B + C the squared center-of-mass energy s = 2EraA + rn^ must be 

(m_B + mcY _ rnA_ 



bigger than (ms + mcY , leading to Eth 
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On the other hand, although the direct experimental upper limits for the effective neutrino 
masses in lepton decays are not very stringent {m^^ < 3 eV, m^^ < 190 keV, m^^ < 18.2 MeV, 
see Ref . |iy] ) , we know that the sum of the masses of light neutrinos (my < niz /2 ~ 45 GeV) 



that have a substantial mixing with z/g, t'^ and Ur is constrained to be smaller than about 



5 eV by their contribution to the total energy density of the Universe . 

The comparison of the cosmological limit on neutrino masses with the energy thresh- 
old in the processes of neutrino detection implies that detectable neutrinos are extremely 
relativistic. This is an important fact that is crucial for the theory of neutrino oscillations. 

Let us consider now the assumption |(A2)| . The flavor state |z/q,) in Eq. (|l]) is defined in 



order to be annihilated only by the flavor field VaL- 

{OWaLWp) OC 5af3 ■ (7) 

Let us check if this is true. The Fourier expansion of the left-handed components of the 
quantized fields of massive neutrinos is (see, for example, |ll2| ) 



ip-x 



where h is the helicity. The neutrino and antineutrino destruction operators a\l^\p), b^j^\p) 
satisfy the canonical anticommutation relations (for Majorana neutrinos af'\p) = b\^\p)) 

{a^k\p) , 4'^\P')} = {bf\p) , bf'^\p')} = 2E {2'nf5\p - p') 6,,, 5,, . (9) 

Hence, 



{OWkL{0)WAP,h))=ufl{p)Sk,, (10) 



and we haveQ 



{oKLmMp^h)) = J2UakU;^{oWkLm^Ap^h)) = Y.u^kU;,ufl{p) ^ (ii) 

k,j k 

However, since detectable neutrinos are extremely relativistic, the contribution of neutrino 
masses in Eq. (|TT]) can be neglected, leading to 

(0|z/aL(0)|z/^) oc 5^f/,fcf/^, = (12) 

k 

because of the unitarity of the mixing matrix. Therefore, the flavor states (^ describe 
correctly neutrinos produced in weak interaction process only in the extreme relativistic 
approximation [113], M, which is valid in neutrino oscillation experiments. Notice that the 



^ For simplicity, we consider a flavor neutrino with helicity h and definite momentum p, according with 
assumption |(A3)| . It is clear that the same conclusion is reached considering different momenta for the 
massive neutrino states. 
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flavor states (|^) are not quanta of tlie flavor fields z/^, only appropriate linear combinations^ 
of the massive neutrino states quanta of the fields u^. 

In order to discuss the physics of neutrino oscillations it is convenient to consider the 
simplest example of neutrino production: pion decay at rest, 

7r+ -> /i+ + z/^ . (13) 

The energy and momentum of each massive neutrino emitted in this process is determined 
by energy- momentum conservation [|T71: 



Since detectable neutrinos are extremely relativistic, only the first order approximation in 
the mass contribution is relevant: 



Pu^E-^^, Eu^E + {l-Oi^, (16) 



with 



Although the relations (|16D have been derived in the specific case of neutrinos produced 
in pion decay at rest, they are valid in general for any process. Indeed, the first order 
approximation in the mass contribution must be proportional to because of the relativistic 
dispersion relation (0). In order to get a quantity of dimension energy, m\ must be divided by 
E, that is the only available energy. The values of E and ^ are determined by the production 
process. 

From Eq. ( |I6D it is clear that in general the equal momentum assumption KA3)| does 
not correspond to reality, unless one considers a special production process in which ^ = 
0. However, as we will see in the following, the oscillation probability turns out to be 
independent of ^. Hence, the incorrect equal momentum assumption = 0) leads to the 
correct transition probability (H). 

How to take into account the different momenta of massive neutrinos in the derivation 
of the oscillation probability? The solution to this question starts from noticing that the 
oscillation probability should be Lorentz invariant because different observers measure the 
same flavor transition probability. But the probability (^ is not Lorentz invariant. In order 
to obtain a Lorentz invariant oscillation probability it is necessary to take into account not 
only the time evolution of the massive neutrino states, given in Eq. (^, which depends on 
their energy, but also their evolution in space, which depends on their momentum. 



See Refs. ^ for a different point of view. 
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The state describes a massive neutrino at the production point x = at the produc- 
tion time t = 0. The state that describes the same massive neutrino at the coordinate x at 
the time t is obtained by acting on lu^) with the space-time translation operator e"*-^''^^': 

|z/fc(x,t)) = e-*^'=*+*^'=nz/fc). (18) 

Equations (|lD and ([18D lead straightforwardly to the Lorentz invariant oscillation probability 



iL,T) 



ETJ* JPkL-iEkT Tj 



13k 



(19) 



which depends on both energy and momentum of massive neutrinos. The probability (|T^ 
describes oscillations in space and time. In order to obtain the probability of oscillations in 
space it is necessary to express the time T in terms of the distance L. 

In connection with the probability ([19|) , it has been claimed recently by some authors 
T9| , pO| ] that the equal momentum assumption |(A3)| should be replaced by the "equal 



energy assumption" 

(A3') The massive neutrino states \ have the same energy = E, but different momenta, 



Pk = ^E^-ml:^E-^. (20) 

This assumption could appear to be attractive because it leads to the vanishing of the time 
dependence of the probability (|1^), which becomes a probability of oscillations in space. 

From Eqs. ( |T6| ) and (p!7| ) one clearly see that in general the equal energy assumption 
is incompatible with energy-momentum conservation in the production process. Therefore, 
unless one considers a special production process in which ^ = 1, the equal energy assumption 
does not correspond to reality. However, since the oscillation probability turns out to be 
independent of ^, the incorrect equal energy assumption = 1) leads to the correct transition 
probability (|^), as well as the incorrect equal momentum assumption (^ = 0). 

There is another simple argument that shows that the equal energy assumption, as well as 
the equal momentum assumption, in general do not correspond to reality: Lorentz invariance 
implies that even if different massive neutrinos have the same energy (momentum) in one 
Lorentz frame, they have different energies (momenta) in all the other frames boosted along 
the neutrino propagation path . 

Indeed, let us assume for example that in a Lorentz frame 5* different massive neutrinos 
have the same energy Ek = E, independent from the mass index k. In this frame the 
momenta of the massive neutrinos are given by Eq. (|20|). 

In another Lorentz frame S" with velocity v with respect to S along the neutrino path 
the energy of the k^^ massive neutrino is 

2 2 

E'k = l{E + vp,)c:.^{l+v)E-^v^ = E'-^^, (21) 
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where 7 = (1 — f ^) and E' = \Jyz^ E is the energy of a massless neutrino in S'. The 
difference between the energies of the /c*^ and j^^ massive neutrinos in the frame S' is 



A£i,^£;-£^ = -^^. (22) 

For relativistic velocities [v ~ 0.1 — 1), the energy difference is of the same order as the 
momentum difference, 

Api^^p',-p'^ = -^^. (23) 

Therefore, it is clear that in the Lorentz frame S' the energies of different massive neutrinos 
are different and the equal energy assumption is untenable. 

Is the transformation from a Lorentz frame S* to a frame S' moving with relativistic 
velocity with respect to S important in practice? The answer is yes. Let us consider, for 
example, the simple case of pion decay (pTSI). For the sake of illustration, let us consider 
the equal energy assumption to be valid for pion decay at rest, even if this assumption 
is incompatible with energy-momentum conservation, as discussed above. Then S is the 
Lorentz frame in which the pion is at rest. 

Many experiments measure the oscillations of neutrinos produced by pion decay in flight. 
These are short and long baseline accelerator experiments and atmospheric neutrino exper- 



iments (see p2[ for a review). The energy of the pions goes from a few hundred MeV (for 



example in the short baseline accelerator experiment LSND ||2^) to hundreds of GeV (for ex- 



ample in the upward-going events measured in the Super-Kamiokande atmospheric neutrino 



experiment ^4 



It is clear that even if the equal energy assumption is valid for pion decay at rest, it cannot 
be valid even approximately in the case of short and long baseline accelerator experiments 
and atmospheric neutrino experiments. Indeed, considering for example a neutrino emitted in 
the forward direction by a pion decaying in flight with energy E.,r ~ 200 MeV, the laboratory 
frame S' is boosted with respect to the frame S in which the pion is at rest by a velocity 
V ~ 0.71, which gives 

— ^ ~ 2.4 , ~ 3.4 . (24) 

1 — V 1 — V 



From Eqs.(^) and ( [23| ) one can see that the energy and momentum difference between 
different massive neutrinos is of the same order of magnitude. Obviously, increasing the pion 
energy, the energy and momentum differences increase and tend to the same limit. 

Let us emphasize that one would obtain the same result choosing another Lorentz frame 
in which the energies of different massive neutrinos are assumed to be equal: from Lorentz 
invariance the equal energy assumption cannot be simultaneously valid for all neutrino os- 
cillation experiments in which neutrinos are produced by pion decay and it cannot be even 
valid in one experiment in which the decaying pion have a spectrum of energies (as always 
happens in practice). 
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Figure |T| 



Figure ^ 



Another obvious problem of the equal energy assumption, as well as the equal momentum 
assumption, is the arbitrariness of the choice of the Lorentz frame in which it is valid, which 
is not based on any physical argument. 

Since the equal energy and equal momentum assumptions are incompatible with Lorentz 
invariance and energy-momentum conservation in the production process, it is better to for- 
get them and consider the different energies and momenta of the massive neutrinos given 
by energy-momentum conservation in the production process. Using the relativistic approx- 
imations in Eq. (p!6D, the probability (|19D becomes 

P,^^,^(L,T) = 5^|f/,fcn?7^,|2 + 2Re5^f/:,f/^,f/«,?7^^.e-^«— (25) 

k k>j 

that describes oscillations in space and time, which depend on the characteristics of the 
production process through the quantity C,- 

In order to obtain the probability of oscillations in space, it is necessary to express 
the time T in terms of the source-detector distance L. This cannot be done considering 
massive neutrinos as plane waves, which extend over all space at all times. It is necessary to 
describe massive neutrinos as wave packets which are localized in the production region at 
the production time, propagate for a distance L during the time T and are localized in the 
detection region at the detection time. 

The production, propagation and detection of a particle wave packet is illustrated by the 
space-time diagram in Fig. |1|. The production and detection processes must be localized, but 
cannot be points in space-time, because in that case the uncertainty principle would imply 
that the energy and momentum have infinite uncertainty. Therefore, the production and 
detection processes occur coherently in finite space-time regions whose sizes are connected 
to the energy-momentum uncertainty by the uncertainty principle. 
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If the production process P occurs coherently in a finite space-time region of dimensions 
(Axp, Atp), the emitted particle is described by a wave packet with size given by 

al ~ Axl + Atl , (26) 

as illustrated in Fig. |I|, in which the space-time coherence region of the detection process 
has dimensions (Axd, Atn)- In Fig- H the propagating particle is assumed to be massless, 
leading to the relation T = L. 

If only one particle propagates between the production and detection processes, as in 
Fig. I, the propagating particle wave packet overlaps with both the production and detec- 
tion processes. In this case the wave packet description does not lead to any interesting 
consequence and can be replaced by a plane wave treatment for all purposes except the 
derivation of the relation between T and L. On the other hand, in neutrino oscillation ex- 
periments the propagating neutrino is a coherent superposition of massive neutrinos, each 
described by a wave packet, as illustrated in Fig. ^ in the simplest case of two massive neu- 
trinos with mi = and m2 > 0. The group velocity of each neutrino wave packet depends 
on its mass: 

Different wave packets are emitted simultaneously by the production process, but arrive at 
the detection process at different times 



ml 



The time T in the phase 



Am?- Am?, 
$..(^, T) = -i^^ - (1 - -i^T (29) 

of the interference term of the fc**^ and j**^ massive neutrinos in the oscillation probability 
(^) must be averaged in the interval 

AT~[t,,4] = L + ^L±^L, (30) 

when the wave packets of v-^ and Vj overlap with the detection process (for m^ > m^). In 
Eq. (^) we have defined Sm^.^- = m^ + m^. The value of the phase (^) in the time interval 
(^) is given by 

, ^ Am?.L ^ Am?.L /Em?. ± AmL\ , ^ 

t.,(L).-^-(l-0^(^^^). (31) 

The second term on the right-hand side of Eq. (|3TD is strongly suppressed with respect 
to the first one because detectable neutrinos are extremely relativistic (Sm^^ <^ E'^ and 
Amlj <^ E"^). Moreover, in real experiments the oscillations due to Am^j are observable only 
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if ~ 1 because of the average over the neutrino energy spectrum, the energy resolution 
of the detector and the source-detector distance uncertainty. Therefore, the leading term 
—Am1jL/2E must be of order one and the correction due to the second term on the right- 
hand side of Eq. (|3lD is negligible. In other words, if oscillations due to Amlj are observable 



in a real experiment, the phase ^kj is practically constant in the time interval [tj , tk\- 

^ 2^ ■ (^2) 

As one can see from Eq. (^, this is the standard value for the phase due to the interference 
of the fc*'^ and j**^ massive neutrinos. It is important here to notice that: 

A. The wave packet description of massive neutrinos is necessary to justify the approxi- 
mation T ~ L (standard assumption |(A4 )| ) that leads to the standard phase (|32D. 



B. The quantity ^ has magically disappeared from the phase, thanks to the relativistic 
approximation and the approximation T ~ L. This is very important because it implies 
that neutrino oscillations do not depend from the specific details of the production 
process. 

As one can see from Fig. 0, another important consequence of the wave packet description 
is that the wave packets of different massive neutrinos propagate with different velocities 
and tend to separate (it is possible to show that the spreading of extremely relativistic wave 
packets is negligible). At a distance larger than the coherence length U^^ the wave packets 
of the fc*'^ and j**^ massive neutrinos cannot both overlap with the detection process and 
the corresponding interference term is suppressed [^. The coherence length L^"*^ can be 



estimated by equating the separation of the wave packets, 

|Am? .| 

lAxfcjl = - T ^ L , (33) 

to the maximal separation allowed for interference, 

I AxlLx -^1 + + A^D ~ ^4 + A4 + ^^l + A^D , (34) 

leading to 



' coh 



I Ax I max ■ (35) 



In conclusion we would like to remark that: 

I. The standard expression (^) for the neutrino oscillation probability in vacuum is robust. 

II. The relativistic approximation is crucial in order to obtain an oscillation probability 
that does not depend on the specific details of the production and detection processes. 

III. The oscillation probability must be Lorentz invariant because different observers mea- 
sure the same flavor transition probability. 
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IV. The equal momentum or energy assumptions do not correspond to reality. They are 
incompatible with Lorentz invariance and with energy-momentum conservation. We 
have shown that they are not needed for the derivation of the oscillation probability 
in space and time. 

V. The wave packet treatment is necessary in order to understand the physics of neutrino 
oscillations and the approximation T ~ L that allows to obtain the measurable oscilla- 
tion probability in space from the Lorentz invariant probability of oscillations in space 
and time. 
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